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Abstract 

Starting from the coherent state representation of the evolution operator with the help of 
the path-integral, we derive a formula for the low-lying levels E = €q — 2Aecos(s + .^)7r of a 
quantum spin system. The quenching of macroscopic quantum coherence is understood as 
the vanishing of cos(s + 0^ disagreement with the suppression of tunneling (i.e. Ae = 0) 
as claimed in the literature. A new configuration called the macroscopic Fermi-particle is 
suggested by the character of its wave function. The tunneling rate (^^) does not vanish, 
not for integer spin s nor for a half-integer value of s, and is calculated explicitly (for the 
position dependent mass) up to the one-loop approximation. 



The question whether or not superpositions of macroscopically distinct quantum states 
exist in nature has attracted considerable attention in recent years owing mainly to the de- 
velopment of technology in mesoscopic physics. Thus Leggett et al.|^ predicted that the 
most intriguing quantum effect which could take place on the macroscopic scale is quantum 
tunneling[Q. Macroscopic magnetisation tunneling is a subject which is being investigated 
extensively and is of growing interest 0, ^]. In discussing this, it is essential to distinguish be- 
tween macroscopic quantum tunneling (MQT) and macroscopic quantum coherence (MQC). 
In the case of MQC, tunneling between neighboring degenerate (perturbation theory) vacua 
is dominated by the instanton configuration with nonzero topological charge and leads to 
a level splitting which (apart from a factor of h ) is seen to be the frequency of oscillation 
between, for example, the easy directions (defined by the vacua) of the ferromagnetic parti- 
cle. The tunneling removes the degeneracy of the ground states, and the true ground state 
is a superposition of the degenerate ground states, for instance, the two easy directions of 
the ferromagnetic particle which are macroscopically distinguishable. For the case of MQT 
the tunneling is dominated by the socalled bounce configuration with zero topological 
charge and leads to the decay of metastable states. The quantum tunneling of magneti- 
zation in small ferromagnetic particles has been well established within the leading order 
approximation of the instanton method for various models 0, using the classical Hamilto- 
nian. Recently an interesting phenomenon known as topological quenching of MQC has been 
studied extensively 0; H, El- The quenching of MQC there was analysed as the suppression 
of tunneling (i.e. the vanishing of the level splitting Ae as given by our eqs.(14,15) below 
which gives its relation to the transition amplitude ) by destructive interference of Euclidean 
paths in the barrier [Q. In the case of MQC the Hamiltonian has a periodic nature, that 
is, there is an infinite number of degenerate vacua by successive 27r extension. However, 
tunneling in the ferromagnetic particle has only been considered in the literature so far as 
that for the double-well potential, like the (p'^ potential in field theory, where the side bar- 
riers are infinitely high. We argue here that the method of calculating the tunneling effect 
for a periodic potential |]10[ ought to be used in dealing with MQC, regarding the periodic 
potential as a one-dimensional super-lattice. The tunneling through one barrier leads to 
the level splitting which extends formally to an energy "band" by translation symmetry (the 
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rotation symmetry in the present case, V{(f) + nvr) = V{(f))). The energy "band" is formally 
the same as that of a one-dimensional tight-binding model in solid state physics. If one 
proceeds like this, an interesting observation emerges. The energy "band" for half-integer 
spin s is quenched to a single degenerate level (doubly degenerate in the present case) in 
agreement with Kramer's theorem which says that for half-integer spin s the degeneracy 
cannot be removed in the absence of the magnetic field. The quenching of MQC is then 
understood as due to the quenching of the energy band, and is the result of the destructive 
interference of Bloch waves. From a theoretical point of view it may be interesting to note 
that the wave function of the giant particle is a superposition of macroscopically distinct 
states with a phase resulting from the Wess-Zumino term. This phase changes its sign un- 
der a 27r-rotation. The ferromagnetic particle thus considered can therefore be regarded 
as a macroscopic Fermi-particle configuration. To see the difference of the two formalisms 
explicitly we start in the following from the spin coherent state path integral. 
The Hamiltonian operator H of the ferromagnetic particle we consider is 

H = KiS,^ + K2Sy (1) 

where Si, i = x,y, z, are spin operators with the usual commutation relation [Si, Sj] = ieijkSk 
(using natural units throughout) and Ki > K2 > 0. This model describes XOF-easy plane 
anisotropy with an easy axis along the x direction We begin with the evaluation of the 
matrix element of the evolution operator in the spin coherent-state representation by means 
of the coherent-state path-integral, i.e. 

< nj\e-'^^^\ni >= J {UtL-,'dfiin,)}{UtU < n^l : e''^^ : \n,^^ >} (2) 

Here we define \nM >= |nj >,|no >= |nj > and tf — ti = 2T, e = ^ respectively. Further 
|n > denotes the spin coherent-state generated from the reference spin eigenstate |s,s > 
such that 

|n >= e^*^— ^^+|s,s >, (3) 
where ( = |e~"^ with polar angle 9 and azimuthal angle of the unit vector n = (sin 9 cos 0, 



sin 6' sin 0, cos 0). The measure of integration is defined by [11 



2s + 1 

dfi{nk) = — duk, diik = sm9kd9kd(f)k (4) 
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The normal ordered infinitesimal evolution operator is 

: exp{—ieH} := exp{ie^ S'^}exp{—ie[-^ S"^ + {Ki ^) Sl]}exp{ie^ S"^} (5) 

The element in eq. (2) can be evaluated by regarding the coherent state formally as an 
"eigenstate" of operators S'^,Sz and 5*4. with "eigenvalues" s{s + l),scos6' and se*'^sin^ 
respectively, the operators in eq. (2) then being replaced by their corresponding "eigenvalues" 
and the coherent-state path-integral in eq.(2) is formally written as the path integral in phase 
space for large s, i.e. 



with canonical variables (p and p = scos6, where 

JZ = ^p-H{<p,p) (7) 
is the phase space (or first order) Lagrangian, and the Hamiltonian 

is just the classical Hamiltonian as in refs. [1^. In H the mass 
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m{(f)) 



2Ki(l - A sin 
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is position dependent (reflecting the fact that there is a velocity dependent force) and 
V{(f)) = K2S^ sin^ (j) is the potential, and A = ^ is a dimensionless parameter. The po- 
sition dependent mass also implies that the original space on which the continuous classical 



mechanics is defined for the spin system (see the following) has an intrinsic curvature flS 
The position dependent kinetic term corresponds to a velocity dependent potential, and this 
requires that one start from the phase space path-integral as pointed out long ago The 



phase of the prefactor e '^^^'^f "^^^ which originated from the inner product of coherent states 
[0,[Tl,i.e. 

< nfc|nfc_i > 

can be put into the Lagrangian, i.e. the angle difference 0/ — 0j is written as 0/ — 0j = J^. (f)dt, 
and has been identified formally as a Wess-Zumino term p, |T^. This phase factor, which 



depends only on the difference of initial and final azimuthal angles of the states |nj > and 
|nj > will be seen in the following to be just the phase difference of wave functions in the 
coherent state representation. We remark here that one cannot add 2?™ (as in eq. (4) of 
ref.Q) into the angle difference since the equality with the original phase factor e~^^^f~'^^^'^ 
would be violated when s is a half-integer. This may be the technical point in ref.0 which 



leads to the vanishing of tunneling. It is well known |]13|, |16| that the appropriate space 
on which the continuous classical mechanics is defined for the spin system is not the unit 
sphere S"^ which is simply connected but the group manifold of S0{3) = SU{2)/Z2 which 
is doubly connected. The path integral quantization of the continuous classical model on 
the group manifold of 5*0(3) gives rise to two inequivalent sectors corresponding to integer 
and half-integer spin. The spin coherent-state representation leads to the path integral 
on S*^. A Wess-Zumino term arises naturally to take care of the topological change of 
space. Integrating out the momentum in the path-integral we obtain the configuration 
space functional integral, i.e. 

< ny|e-2*^^|n, >= e-'^^f-'^^>]C{<j)f,tf = T;<P„U = -T) (9) 

where 

}Ci<P,,tr,M= fv^Pe^l''^''^''^'' (10) 



The functional JC is the Feynman propagator in configuration space with the second order 
Lagrangian 



and the measure 



P0 = n^i-/i 



lm{(f)k 



27[ie 

The canonical momentum is defined as usual, i.e. p = m{(f))(j). The dynamical system is 
then reduced to that on the circle 5*^. The potential is periodic with period vr and has an 
infinite number of degenerate minima by successive 2tt extension at $„ = ravr, n being an 
integer. Considering the periodic potential as a super-lattice with lattice constant vr, we can 
derive the energy spectrum in the tight-binding approximation. We let |m, $„ > be the set 
of low-lying eigenstates of the zero order (i.e. oscillator-approximated) Hamiltonian in 
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the nth well near the local minimum such that 

Ho\m,<^„.>=^m\m,<^n> (11) 

and we define lip > as the eigenstate of the full Hamiltonian with translation symmetry, 
H{(j)) = H{(f) + nvr) and H\%1) >= E\iIj >. The state {ip > is then the (unnormalized) Bloch 
state 

|V>=S„e*«*"|m,$„> (12) 

Here {ip > and |m, $„ > denote the usual Dirac kets. Eqs. (11) and (12) just define the 
ordinary quantum mechanical problem where ^ is the Bloch wave number. Since the range 
of the azimuthal angle is 27r, or in other words the physical length of the superlattice is 
27T, the periodic boundary condition 

implies that ^ must be an integer. The character of spin comes into the picture only through 
the spin coherent-state representation which leads to a topological phase of wave functions 
resulting directly from the Wess-Zumino term. The energy "band" is seen to be given by 

E = - 2Aecos(s + 07r (13) 

where the expression Ae is given by 

Ae = - y <(0, <^n)HUm{<P, $n+l)# (14) 

which is the usual overlap integral or simply the level shift due to tunneling through any 
one of the barriers where Mm(0, $n) denotes an eigenfunction of the Hamiltonian Hq for 
the nth well. The energy formula eq. (13) is only formally the same as that of the energy 
"band" of a one-dimensional lattice. In our case there are only two levels for integer s and 
there is only one level for half-integer s. The instanton method of calculation is nothing 
but a procedure to evaluate the overlap integral analytically. If s is half-integral the energy 
"band" of eq.(13) shrinks to a single level which is doubly degenerate [|l^ - in agreement 
with Kramer's theorem and with numerical calculations 0. 
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However, the quenching of MQC is due to the vanishing of cos(s + 0'"" to the 

vanishing of tunnehng (i.e. Ae = 0) as conjectured in ref.|P]. 
Since ranges from to 27r, the Bloch waves 

(remembering that in the spin coherent-state representation the state lip > becomes the wave 
function < (plip >= ip[(f))e~^'^'^) then obey the periodic boundary condition for integer spin, 
and the antiperiodic boundary condition for half-integer spin. Therefore the wave function 
of the giant spin is a superposition of macroscopically distinct states and changes sign under 
rotation of 2tt for a half-integer spin s. The ferromagnetic particle with half-integer spin s 
may be looked at as a macroscopic Fermi particle. Up to this point we worked only in real 
time. The above results do not depend on Euclidean paths inside the barrier in disagreement 
with the formalism in ref.0. 

We now calculate the tunneling parameter, i.e. the level splitting Ae and hence the 
tunneling rate which does not vanish for half-integer spin s. Passing to imaginary time 
by Wick rotation t = it, (3 = iT, the amplitude for tunneling from the initial well, say that 
with n = (and (pi = 0), to the neighboring well with n = 1 (and 0/ = vr) and considering 
large (3, the transition amplitude between the corresponding coherent states is seen to be 

< n(7r)|e-2/5-f^|n(0) > = < n(7r)|0, $i >< 0, $o|n(0) > e'^^^" sinh(2/?Ae) 

(15) 

where 

is the Euclidean propagator with Euclidean action defined by 

Se = f CEdr, Ce = \m{ct>)^'' + V{<t>) (16) 

In eq.(15) cq is the ground state energy of Hq. From now on = ^ denotes the imaginary 
time derivative. The phase in eqs. (6) and (9) for the case 0/ = vr, 0j = is just the phase 
relating to wave functions in the 0th and 1st wells in the coherent state representation. 
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In the following the Euclidean propagator /C^; is evaluated with the instanton method. 
After evaluation we compare the result with eq.(15) to find the level splitting Ae. The 
instanton configuration which minimizes the Euclidean action Se is 



(/)c = arcsin[cosh ti;o(r — tq) — A sinh ujo{t — tq)] 2 (17) 

with position tq and cJq = AKiK2s'^. The Euclidean action evaluated for the instanton 
trajectory eq.(17), sometimes called the instanton mass, is 

/CO . 1 _|_ ^/T 

m{<P,)<f)ldT = s In ---^ (18) 
-00 1 — V A 

in agreement with refs.|^, [1^. The functional integral ICe can be evaluated with the station- 
ary phase method by expanding about the instanton trajectory 0c such that (p = 4>c + r], 
where r] is the small fluctuation with boundary conditions ?](/?) = rj{—[3) = 0. Up to the 
one-loop approximation we have 

ICe = e^^^I (19) 



where 



/ = / Vr]e-^^^ (20) 

Jvi-f^)=o 



is the fluctuation functional with the fluctuation action 

SSe = f vMvdr (21) 

where 

M = -^^m(0,)^ + ^(0j (22) 

with 

\/(0,) = l[-m'(0,)0, - ^m"(0,)0^ + V"{<P,)] (23) 

Here M{(j)c) is the operator of the second variation of the action and m'{(j)c) = ^^^^ I ^=ci>c ■ 
As in the usual method of evaluating the fluctuation integral /, we expand the fluctuation 
variable rj in terms of the eigenmodes of M and set rj = Y^nCnipni where ipn denotes the nth 
eigenfunction of M, and express the result of the integration as an inverse square root of 
the determinant of M. In view of the time translation symmetry of the equation of motion, 
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the functional integral /C^; is not well defined when expanded about the classical solutions. 
The translational symmetry results in zero eigenmodes of the second variation operator M 
of the action (which in the present case, of course, has only one). This problem can be 



cured by the Faddeev-Popov procedure |T8[ or in a more systematic way with the help of 
the BRST transformation |]19|. Then the zero-mode functional integral is converted into an 
integral over the collective coordinate, i.e. the instanton position tq, and leads to 2(3 and a 
Faddeev-Popov determinant ^/D with 

<-) 

The complicated determinant of M and the related Jacobian are avoided by the direct 
integration involved in the transformation of the well-known shift method [|T9l , i.e. 

v = y + kr ^M^W (25) 

Jn 02 (r') 

where the derivative of the instanton solution, i.e. 0c, is simply the zero eigenmode of M. 
The resulting fluctuation integral is 

^ ^ [kmd-f^) r . .J ~^ (26) 



V27r J-fi (f)l{T)m{(j)c{T)) 

Following the procedure of ref.|T9| or earlier work of ref.[^] the one instanton contribution 
to the propagator in the one-loop approximation is calculated to be 

^« ^ 2(3^(1 - A)-^s2i^2e-^o^e-^^ (27) 

71 

To obtain the desired result proportional to sinh(2/5Ae) in eq.(15), the contributions of the 
infinite number of instanton and antiinstanton pairs to the one instanton contribution have 
to be taken into account. Interactions among instantons and antiinstantons are neglected 
in the dilute instanton-gas approximation. The contribution of one instanton plus n such 
pairs to the propagator is obtained with the help of the group property of the Feynman 
path-integral, i.e. 



Summing over all contributions /C^"^^\ the final result of the propagator is found to be 

]Ce = A3(^)5e-/^-o sinh[2/3.2t{^^^^}h^sie-^=] (29) 
27r (1 — AjTT 

Comparing with eq.(15) the level splitting is seen to be 

Ae = 22{-- —^2XAs-2e (30) 

(1 - A)7r 

The product of the values of the wave functions at neighboring potential minima is 

< n(7r), 0|0, $1 >< 0, $o|n(0) >= {^)hh-'''' (31) 

and the ground state energy of Hq is identified as eo = ^ujq- Expanding Ae in the region of 
small values of A, the level splitting is 

Ae = 25{ — L^}2A4s2e-2^^^ (32) 

TT 

Since s is a large number (500 to 5000 as cited in the literature 0]), the level splitting is 
suppressed in leading order by the factor 

For Ae to be observable, plausible values of A are < 10~^ as discussed e.g. in ref.0. 
The quantum correction proportional to ei''^'^ cannot increase the tunneling probability 
significantly, but may be of importance for small values of A. We remark that one of the 
crucial points of the effective Hamiltonian as well as the Lagrangian for a quantum spin 
system is the position dependent mass which requires great care in the evaluation of the 
path integral. The result is not simply the factor y/B^ as in the usual case of constant mass. 

If in the case of the potential with degenerate minima only one of the potential minima is 
considered as a perturbative vacuum; however, a fictitious imaginary energy can be calculated 
by consideration of possible back and forth tunneling in the barrier. Therefore there exists 
a relation between the level splitting and this imaginary energy, and has been referred to as 



the Bogomolny-Fateyev relation ||T0|, P3| 



ImE = 27T(Ae? 
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in the context of large order perturbation theory (in [|T^, ^ the factor 27r arises with 
ujq = 1). Recently the relation has been rederived and generalized to low-lying levels with 
the instanton method by regarding the instanton-antiinstanton pair as a bounce-like con- 
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In the literature the quenching of the energy "band" structure for half-integer s (here, 
of course, we have only formally a "band") has been interpreted as due to the quenching of 
tunneling. The level splitting or overlap integral is said to vanish due to an interference of the 
Euclidean paths with the topological phase in eqs.(6) and (9). The destructive interference, 
however, does not occur for closed loops because the interference in this case would be 
proportional to |e'^'^'* + e~*^'^*| which does not vanish for half-integer values of s. Thus 
there is an apparent contradiction with the Bogomolny-Fateyev relation from the viewpoint 
of quenching of tunneling: the level splitting is said to vanish but its theoretical partner, 
namely the imaginary part which is calculated from the bounce (a closed Euclidean path) 
does not. Therefore we suspect a misinterpretation in ref. as we indicated in our earlier 
remark. There is no contradiction in our framework. The level splitting or overlap integral 
does not vanish, which is, as it ought to be. 

In conclusion we can say: In a theoretical analysis of MQC of ferromagnetic particles, 
the periodic nature of the Hamiltonian must be taken into account. The quenching of MQC 
results from interference of Bloch waves which obey an antiperiodic boundary condition. 
The tunneling cannot be suppressed for half-integer spin s in a consistent instanton method. 
The ferromagnetic particle with half-integer spin s may be looked at as a macroscopic Fermi 
particle configuration. The correction of quantum fluctuations with the position dependent 
mass to the level splitting for MQC, here first reported ||2^, does not enhance the tunneling 
probability significantly, as expected. Nevertheless the contribution of quantum fluctuations 
increases the tunneling probability along with s, which is proportional to the volume of the 
giant particle. In the above eqs. (30) and (32) provide explicit formulas of the level splitting 
suitable for quantitative evaluation. Enhancement of MQC may take place theoretically if 
tunneling from an excited state is considered, even for low-lying levels. This investigation 
will be reported elsewhere. 
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